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O \ Abstract 

, _. ' Continuous limits of discrete systems with long-range interactions are considered. 

The map of discrete models into continuous medium models is defined. A wide class 
of long-range interactions that give the fractional equations in the continuous limit is 

i-Q . discussed. The one-dimensional systems of coupled oscillators for this type of long- 

I , range interactions are considered. The discrete equations of motion are mapped 

'^ \ into the continuum equation with the Riesz fractional derivative. 



^ '■ 1 Introduction 



(N 



Equations which involve derivatives or integrals of noninteger order [H [21 [3l HI [5] have 
found many applications in recent studies in mechanics and physics [SI El El El [IDl E] • 
Usually the fractional equations for dynamics or kinetics appear as some phenomenological 
Q \ models. Recently, the method to obtain fractional analogues of equations of motion was 

considered for sets of coupled particles with a long-range interaction [121 [131 [H] . Examples 
k> \ of systems with interacting oscillators, spins, or waves are used for numerous applications 

^ ■ in physics, chemistry, biology [I51IISllI7llISllIHll2ni[2Il[22l[2ni[2a[2Sl[2S]. Transfer 

from the equations of motion for discrete systems to the continuous medium equation 
with fractional derivatives is an approximate procedure. Different applications of the 
procedure have already been used to derive fractional sine-Gordon and fractional wave 
Hilbert equation [121 HI], to study synchronization of coupled oscillators [13j, and for 
fractional Ginzburg-Landau equation |13j . 

Long-range interaction has been the subject of great interest for a long time. Ther- 
modynamics of the model of classical spins with long-range interactions has been studied 
in Refs. [T5l [T6l [T9l [T7] . An infinite one-dimensional Ising model with long-range in- 
teractions was considered by Dyson [ISJ. The d-dimensional classical Heisenberg model 



with long-range interaction is described in Refs. ^E\ ^U\ , and their quantum generahza- 
tion can be found in Ref. [17]. The long-range interactions have been widely studied 
in discrete systems on lattices as well as in their continuous analogues. Solitons in a 
one-dimensional lattice with the long-range Lennard- Jones-type interaction were consid- 
ered in Ref. [27]. Kinks in the Frenkel-Kontorova model with long-range interparticle 
interactions were studied in Ref. [2H]- The properties of time periodic spatially localized 
solutions (breathers) on discrete chains in the presence of algebraically decaying interac- 
tions were considered in Refs. [211 [25|. Energy and decay properties of discrete breathers 
in systems with long-range interactions have also been studied in the framework of the 
Klein-Gordon [22|, and discrete nonlinear Schrodinger equations (22]. A remarkable prop- 
erty of the dynamics described by the equation with fractional space derivatives is that 
the solutions have power-like tails. Similar features were observed in the lattice models 
with power-like long-range interactions [23] [211 1251 [301 1311 ISll |ll|. As it was shown in 
Ref. [131 El, analysis of the equations with fractional derivatives can provide results for 
the space asymptotics of their solutions. 

The goal of this paper is to study a connection between the dynamics of system of 
particles with long-range interactions and the fractional continuous medium equations by 
using the transform operation. Here, we consider the one- dimensional lattice of coupled 
nonlinear oscillators. We make the transform to the continuous limit and derive the 
fractional equation which describes the dynamics of the oscillatory medium. We show 
how the continuous limit for the systems of oscillators with long-range interaction can be 
described by the corresponding fractional equation. 

In Sec. 2, the equations of motion for the system of oscillators with long-range interac- 
tion are considered. In Sec. 3, the transform operation that maps the discrete equations 
into continuous medium equation is defined. In Sec. 4, the Fourier series transform of 
the equations of a system with long-range interaction is realized. In Sec. 5, we consider a 
wide class of long-range interactions that can give the fractional equations in the contin- 
uous limit. In Sec. 6, the simple example of nearest-neighbor interaction is considered to 
demonstrate the application of the transform operation to the well-known case. In Sec. 7, 
the power-law long-range interactions with positive integer powers are considered. In Sec. 
8, the power-law long-range interactions with noninteger powers and the correspondent 
continuous medium equations are discussed. In Sec. 9, the nonlinear long-range inter- 
actions for the discrete systems are used to derive the Burgers, Korteweg-de Vries and 
Boussinesq equations and their fractional generalizations in the continuous limit. In Sec. 
10, the fractional equations are obtained from the dispersion law for three-dimensional 
discrete system. The conclusion is given in Sec. 11. 



2 Equations of motion for interacting oscillators 

Consider a one-dimensional system of interacting oscillators that are described by the 
equations of motion, 

^ = ^J4n) + FK), (1) 

where m„ are displacements from the equilibrium. The terms F{un) characterize an in- 
teraction of the oscillators with the external on-site force. The term Xn{u) is defined 

by 

in{u)= ^ J{n,m) E{Un,Um), (2) 

m=— oo 

and it takes into account the interaction of the oscillators in the system. 
Examples. 

1) If J{n, m) = 6n+l,m - 5n,m, and E{Un, Urn) = Um, then X„(m) = Un+1 -Un = AUn- 

2) For J{n, m) = 6n+i,m - 26n,m + Sn-i,m, and E{un, Um) = Um, we get 

Iniu) = Un+1 - 2m„ + Un~l = ^^Un- 

3) We can consider the long-range interaction that is given by J{n) = \n\~^^'^"\ where a 
is a positive real number. In this case, we have nonlocal coupling given by the power-law 
function. Constant a is a physical relevant parameter. Some integer values of a correspond 
to the well-known physical situations: Coulomb potential corresponds to a = 0, dipole- 
dipole interaction corresponds to a = 2, and the limit a ^ oo is for the case of nearest- 
neighbor interaction. 

For the term ([2]) with E{un,Um) = Um, the translation invariance condition is 

+ 00 

E J(n,m)=0 (3) 

m=— oo 

for all n. If (^ cannot be satisfied, we must define E{un,Um) = Un — Um, and the 
interaction term ([2]) is 

+00 

iniu) = ^ J{n,m) [Un~Um]- (4) 

m=—oo 

This interaction term is translation invariant. Note that the noninvariant terms lead to 
the divergences in the continuous limit (see Appendix). 

In this paper, we consider the wide class of interactions (jl]) that create a possibility 
to present the continuous medium equations with fractional derivatives. We also discuss 
the term ([2]) with E{un,Um) = fiun) — fium) as nonlinear long-range interaction. As the 



examples, we consider f{u) = v? and f{u) = u — gv? that give the Burgers, Korteweg- 
de Vries and Boussinesq equations and their fractional generalizations in the continuous 
limit. 

3 Transform operation 

In this section, we define the operation that transforms the system of equations for Un{t) 
into continuous medium equation for u{x,t). 

To derive a continuous medium equation, we suppose that M„(t) are Fourier coefficients 
of some function u{k,t). We define the field u{k,t) on [—Kl2,Kl2] as 

+00 

U{k,t)= Y, Un{t)e-''^-=J^A{Unm, (5) 

n=— 00 

where Xn = nAx, Ax = 2tt/K is distance between oscillators, and 

1 r+K/2 
Unit) = - dk u{k, t) e*'^^" = r^\u{k, t)}. (6) 

^ J-K/2 

These equations are the basis for the Fourier transform, which is obtained by transforming 
from discrete variable to a continuous one in the limit Ax -^ {K -^ oo). The Fourier 
transform can be derived from ([5]), (El) in the limit as Ax -^ 0. Replace the discrete Un{t) 
with continuous u{x,t) while letting x„ = nAx = 27m/ K -^ x. Then change the sum to 
an integral, and Eqs. ([5]) and ([6]) become 

/ + 00 
dx e-^^^M(x, t) = T{u{x, t)}, (7) 

'OO 

1 f + OO 

u{x, t) = — dk e^^^M(fc, t) = J^-\u{k, t)}, (8) 

where 

u{k,t) = Cu{k,t), (9) 

and £ denotes the passage to the limit Ax — > {K — > oo). Note that u{k,t) is a Fourier 
transform of the field u{x,t), and u{k,t) is a Fourier series transform of Un{t). The 
function u{k,t) can be derived from u{k,t) in the limit Ax -^ 0. 

The procedure of the replacement of a discrete model by the continuous one is defined 
by the transform operation. 

Definition 1. Transform operation T is a combination T = T'^L JF^ oj the opera- 
tions: 
1) The Fourier series transform: 

^A : Unit) -^ J'AiUnit)} = u{k, t). (10) 



2) The passage to the limit Ax -^ 0; 

C: u{k,t) ^ C{u{k,t)} = u{k,t). (11) 

3) The inverse Fourier transform: 

T-^ ■ u{k, t) -^ J^-^{u{k, t)} = u{x, t). (12) 

The operation T = T^^C T/^ is called a transform operation, since it performs a 
transform of a discrete model of coupled oscillators into the continuous medium model. 

Proposition 1. The transform operation T maps the function F{un) into the function 
F{u{x,t)), i.e., 

f F{un{t)) = F{u{x,t)), (13) 

where u{x,t) = Tunit), if the function F satisfies CF{un) = F{Cun). 

Proof. The Fourier series transform leads to 

^A : F{un) -^ J'aFM. (14) 

Note that J^/\F{un) 7^ F{J^^Un) = F{u{k,t)). The passage to the limit Ax — > gives 

C : J^AF{un) -^ CJ^aFM. (15) 

Then 

CJ^A{F{un)} = J'iCFiun)} = ^{F(£m„)} = J^ {F {u{x , t))} , (16) 

where we use CJ^a = J^ C.- The inverse Fourier transform get 

T-^ ■ T{F{u{x, t))} ^ T-^{T{F{u{x, t))}} = F{u{x, t)). (17) 

As the result, we prove flT^ . 

4 Equations for momentum space 

Let us consider a system of infinite numbers of oscillators with interparticle interaction 
that is described by (jl]). We suppose that J{n,m) satisfies the condition 



J{n,m) = J{n — m), y^ \J{n)\^ < 00, 



n=l 



where J{—n) = J{n). 



Proposition 2. The Fourier series transform JF^ maps the equations of motion 
d^unit) 



g ^ J{n,m)[un-Um]+ F{un), (19) 



dt^ 

m=—oo 

where m„, is the position of the nth oscillator, and F is an external on-site force, into the 
equation 

^^^ = glUO) - UkAx)]u{k,t)+J^^{F{ur.)}, (20) 

where 

u{k,t) = J^A{un{t)}, Ja{kAx) = J^A{J{n)}, 

and J-'A{F{Un)} is an operator notation for the Fourier series transform of F{un). 

Proof. To derive the equation for tlie field u{k, t), we multiply Eq. (IT9l) by exp{—iknAx), 
and summing over n from — cxo to +cx3. Then 






^2 +00 +00 +CXD 



n=—oo 






n=—oo m=-oo 


From 






u{k,t) = 


+ 00 

\ ^ -jfcnAx, 
n=—oo 


the left-hand side of 


m 


gives 






+ 00 

n=—oo 


-iknAxd^Mt) 
dt^ 


pa +°o 

dt^ ^ 

n=—oo 



n=—oo 



(21) 

Unit), (22) 

"A.^„(t) = ^!^. (23) 

The second term of the right-hand side of fl2Tl) is 

+ 00 

E e~*'="^^'F(M„) = .FA{F(n„)}. (24) 

n=— oo 

The first term of the right-hand side of fl?Il) is defined by the function J{n,m). Let 
us introduce the notation 

+ 00 

JaikAx) = Y, e-*'="^V(n). (25) 

n=— oo 

Using J{—n) = J{n), the function 1^^ can be presented by 

+ 00 +00 

4 (A; Ax) = ^ J{n) [e^'^'^^^^ + e^^^"^^) = 2 ^ J{n) cos (A; Ax) . (26) 

n=l n=l 

6 



From (12^ it follows that 

Ja{kAx + 27rm) = Ja{kAx), (27) 

where m is an integer. 

The interaction term in (I2T!) is 

+00 +cx> 



Yl Yl e-*'="^V(n,m)K-M„] 



n=— 00 rra=— 00 



hoo +00 +00 +00 



n=— 00 m=— 00 n=— 00 7n=— oo 

Using (I22I) and (!25|) . the first term on the r.h.s. of fl28l) gives 



hoo +00 +00 +00 



J- ^ e-*"^V(n,m)M„= ^ e'^^^^^M^ ^ J(m') = m(A;, t) J„(0), (29) 

n=— 00 m,=— 00 n=— 00 m'=— oo 

where we use f lTS]) and J(m' + n, ra) = J{m'), and 

+00 00 

J^(0) = J] J(n) = 2^J(n). (30) 



n=— 00 ji=l 



For the second term on the r.h.s. of 



hoo +00 +00 +00 



Y Y e"*''"^V(n, m)M^ = ^ ^^ Y^ e-^'^"^V(n, m) 

n=— 00 ?ri=— 00 m=— 00 n=— oo 



-00 +00 



Y Ume-'^'^^'' Y e-"'''''^^J{n) =u{k,t)Ja{kAx), (31) 



m=— 00 n'=— 00 

nVO 

where we use J{m, n' + m) = J{n'). 
As a result, Eq. fl2Tl) has the form 

^^^^ = glUO) - UkAx)]u{k,t)+J'^{F{u^)}, (32) 

where J-'A{F{un)} is an operator notation for the Fourier series transform of F{un)- 



5 Alpha-interaction 

Let us consider the interaction term 



+ CX) 



in{u)= ^ J{n,m)[un-Um]- (33) 

m=—QO 

where 

oo 

J{n,m) = J{n — m) = J{m — n), y^\J{n)\'^ < oo. (34) 

n=l 

In Sec. 4., we prove that the Fourier series transform J^a of fl33l) gives 

J'Aiiniu)} = g[UO) - UkAx)]uik,t), (35) 

where u{k,t) = J-'Aiunit)} , and 

+00 oo 

J^{k)= ^ e-^'='^J(n) = 2^J(n)cos(A;n). (36) 

n=— oo n=l 

Definition 2. T/ie interaction term ^3^ in equation of motion (Q]) zs called a-interaction 
if the function ^3^ satisfies the condition 

,.„^ |J„W - 4(0)] ^ (37) 

fc^o \k\°' 

where a > and < \Aa\ < oo. 

If the function Ja{k) is given, then J{n) can be defined by 

1 r - 

J{n) = — Ja{k)cos{nk) dk. (38) 



TT JO 

The condition fl37j) means that Jq,(/c) — Ja(0) = 0{\k\°'), i.e. 

Ja(A;)- Ja(0) = A,|A;|" + i?,(A;), (39) 

for /c — i> 0, where 

hm /?«(A;)/|fc|° = 0. (40) 

Examples. 

1) The first example of the a-interaction is 



Using (1351) . we obtain 

/(_l)n »+l C_^Nn 1/2 \ 

where Li{ii, i>, z) is the Lommel function p[5j. 
2) The second example of the a-interaction is 

J{n) = L^. (42) 

Using (Ref. [46], Sec. 5.4.2.12) 

E ^ cos(nA:) = i T/c^ - ^ y |/.| < vr, 

71=1 ^ ^ 

we get 

4(A;) = 2 V^^cos(A;n) = -P - -, |A;| < vr. (43) 

Then we have a = 2, and 



n=l 



Ja{k) - Ja{0) = (1/2) k\ (44) 

The inverse Fourier transform of this expression gives the coordinate derivatives of second 
order 

j'-'iUk) - Jam 



I d^ 



2 9x2 

3) For the interaction f H2l) . we have a = 2 and the inverse Fourier transform of (jSj) 
gives the second-order derivative. At the same time, the interaction 



An) 



m? 



gives a = 1 and then the first-order coordinate derivative. It can be proved by using (Ref. 
Hg, Sec. 5.4.2.12) 

L{k) = 2 y ?^^!M. = 1[3A;2 _ evrfc + 2n\ {0 < k < 27r), (45) 

n=l 

and Ja{k) — Ja(0) ~ — ttA; for A; — > 0. Therefore, the inverse Fourier transform leads to 
the derivative of first order. 

4) For noninteger and odd numbers s, 

J{n) = |n|-("+^), s > (46) 



is an a-interaction. 

For < s < 2 (s 7^ 1), we get 

Ja{k) - Ja{0) = 2r(-s) cos(7rs/2) |A;|^ (47) 

For s = 1, 

Mk) - Ja{0) = -(vr/2) k. (48) 

For noninteger s > 2, 

Ja{k)-Ja{0) = ~C{a-l)k\ (49) 

where ({z) is the Riemann zeta-function. The interaction fllBl) is considered in Sec 7. 

5) The other example is 

(-1)" 

J(n) = — ^—4 -. (50) 

^ ' r(l + a/2 + n)r(l + a/2-n) ^ ' 

Using the series (Ref. [16], Sec.5.4.8.12) 

V^ (-1)" . ,x 22/3-1 /^x 1 

^ r(/3 + i + n)r(/3 + i-n) ^°^("') = rWTT) ^^" U J ' 'mP^y ^''^ 

where P > —1/2 and < A; < 2n, we get 

Uk) - i.(0) = fj^^^ sin" Q^ . (52) 

In the hmit fc — > 0, we obtain 

Uk)-U0)^—^—-\kr. (53) 

i (a + 1) 

For noninteger a, the inverse Fourier transform of fl53l) gives the fractional Riesz deriva- 
tive [1] of order a. 



6) The a-interaction 



J(n) 



9 9 ' 



gives 



For A; — > 0, we obtain 



J{k) = /. . cos(aA;) - \. (54) 



asmivra) a^ 



an , 2 



Ja(A:)-Ja(0)- ^ . , . k\ (55) 

2sm aTT 
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The inverse Fourier transform of (155]1 leads to the coordinate derivative of second order. 
7) For J(n) = l/n\, we use 

^^^^M = e^°^'=cos(sinA;), |A;| < oo. (56) 

n=l 

The passage to the hmit A; ^ gives 

Ja{k)-Jc,{<d)^-Aek. (57) 

Then a = 1, and we get the derivative of first order. 

Proposition 3. The transform operation T maps the discrete equations of motion 
d'^Un 



dt^ 

■m=— CO 



g J^ J{n,m) [Un-Urn] + F{Un) (5^ 



with noninteger a-interaction into the fractional continuous medium equation: 

Q2 Qa 

--^«(x, t) - G^A^——u{x, t)~F (m(x, t)) = 0, (59) 

where 9"/9|x|" is the Riesz fractional derivative, and 

G', = ^|Axr (60) 

is a finite parameter. 



Proof. The Fourier series transform JF^ of (!58l) gives (1201) . We will be interested in 
the limit Ax -^ 0. Then Eq. (1201) can be written as 

— m(A;, t) - Go. faAik) u{k, t) - J^a{F (M„(t))} = 0, (61) 

where we use finite parameter (!60|) . and 

faAik) = -^alA^r - Ra{kAx)\Ax\-''. (62) 

Note that R^ satisfies the condition 

i?„(A;Ax) ^ 
lim — — — ; = 0. 

Ax-^O I Ax I" 

The expression for 7^,a(^) can be considered as a Fourier transform of the operator (jlj). 
Note that g —>■ oo ioi the limit Ax -^ 0, if Ga is a finite parameter. 

11 



In the limit Ax -^ 0, Eq. (1^ gets 

—uik, t) - Ga T^{k) u{k, t) - J^{F {u{x, t))} = 0, (63) 

where 

u{k,t) = Cu{k,t), %{k) = Cfo,,A{k) = -v4„|fc|". 

The inverse Fourier transform of (p5]) gives 

--^m(x, t) - Ga Ta{x) u{x, t) - F {u{x, t)) = 0, (64) 

where Ta{x) is an operator 

r^{x)=J^-'{Uk)} = A^^. (65) 

Here, we have used the connection between the Riesz fractional derivative and its Fourier 
transform [1]: 

\k\'"< — ^-T^. (66) 

The properties of the Riesz derivative can be found in Refs. [H E], [3l H] . Note that the 
Riesz derivative could be represented as 

^^(^^ i) = -7, ^—TTT. (^X^, t) + V-_U{^, t)) , (67) 

o\x\°' 2cos(7ra/2) 

where a 7^ 0, 1,3,5..., and V^ are Riemann-Liouville left and right fractional derivatives 
defined by P [21 [3l H 



r(m - a) dx"" J_^ {x - ^)°-'»+i 

where m — \ < a < m. 

As the result, we obtain continuous medium equations (l59l) from ( l64l) and (l65ll . 



6 Simple example of nearest-neighbor interaction 

In this section, we demonstrate the application of transform operation to the well-known 
case: 

J{n, m) = 6n+l,m - 2(5„,„ + 6n-.l,m, (69) 



12 



where Sn,m is the Kronecker symboL Then the interaction term ([2]) has the form 

in{u) = {Un+1 - Un) - (Wn " Un-l) , (70) 

and describes the nearest-neighbor interaction. As the result, equations of motion (1191) 

have the form 

d u 

-—-^ = g[Un+l-2Un + Un-l]+ F{Un). (71) 

The well-known result is the following. 



Proposition 4. The transform operation T maps the equation of motion (71) into 
the continuous medium equation 

^^ = G2^-(-'^) + ^W' (72) 

where 

G2 = g{^xf (73) 

is a finite parameter. 

Proof. To derive the equation for the field ^(fc, t), we multiply Eq. (17T!) by exp(— ifcnAx), 
and summing over n from — oo to +oo. Then 

n=— oo 71=— oo n=— oo 

The first term on the r.h.s. of (|74l) is 

+00 4-00 



n=— oo n=— oo 7i=— oo 

Hoo +00 

— -ji./'cu-i ^ A--, 

m,'=— oo s=— oo 

+ 00 +00 



J2 e-''=('"-i)^"u™-2M(/t,t)+ ^ e~ife(«+i)A. 



iksAx 

etc 



772' = — OO S = — OO 

e"'^''u{k, t) - 2u{k, t) + e-'^^^'uik, t) = [e''^^" + e-"'^'' - 2]u{k, t) = 

= 2[cos (kAx) - l]u{k, t) = -4 sin^ (kAx) u{k, t). (75) 

13 



As the result, we obtain 



g J^{kAx) u{k, t) + T^{F K(t))}, (76) 



where 

JaikAx) = -4 sin^ (kAx) . (77) 

For Ax -^ 0, the asymptotics of the sine is 



°° /_i\m+l 1 



^(2m+l)!~ "~ 6' 



m=0 

and (1771) can be presented by 

Ja{kAx) ^ - {kAxf + -^ (fcAx)^ . (78) 

Using the finite parameter (173|) . the transition to the limit Ax ^ in Eq. (1761) gives 

^^^ = -G2fc^«(fc, t) + ^-^i^M}, (79) 

where we use < IG2I < 00. As the result, the inverse Fourier transform of (1791) leads to 
the continuous medium equation ( 1721) . 

7 Integer power-law interaction 

Let us consider the power-law interaction (jlj) with 

J(n) = \n\-^'+^^ (80) 

with positive integer number s. 



Proposition 5. The power-law interaction (E^ for the odd number s is a-interaction 
with a = 1 for s = 1, and a = 2 for s = 3, 5, 7.... For even numbers s, (E^ is not a- 
interaction. For odd number s, the transform operation T maps the equations of motion 
with the interaction (E3j into the continuous medium equation with derivatives of first 
order for s = 1, and the second order for other odd s. 

Proof. From (l20l) . we get the equation for u{k,t) in the form 

^^^1^ + g lUkAx) - J„(0)] u{k, t) - T^{F {un(t))} = 0, (81) 
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where 

+00 



J«(A;Ax) = Y, e-'^"^"|n|-(^+^). (82) 

n=— 00 

The function fl82l) can be presented by 

+CO ^ +00 ^ 

Ja{kAx) = J2 ^TI (e"'^"^"" + e^^"^^) = 2 ^ ^^ cos (knAx) . (83) 



n=l n=l 



Then we can use (Ref. |16] Sec. 5.4.2.12 and Sec. 5.4.2.7) the relations 

y^-^^ = ^{3e-6nk + 2n'), (0 < A: < 27r), (84) 

n=l 

^^ n2™ 2 2m ! V2vr / ' ^ - - ;, v ; 

n=l ^ ' ^ ' 

where ?ti = 1,2,3,..., and B2m{z) are the Bernulh polynomials [44j. These polynomials 
are defined by 

n 

B^{k) = Y,C'^B,k^-% (86) 

where 5s are the Bernoulli numbers from 

For example, 

B2{k) = k'^-k + l/Q, B4{k) = k^-2k^ + P-l/30. (88) 

Note B2m-i = for m = 2, 3, 4... gl]. 
For s = 1, we have 

Ja{kAx) - Ja(0) = -{kAxf - TTkAx ^ -TTkAx. (89) 

For s = 2m — 1 (m = 2, 3, ...), we have 

'_i \m-l / L 



Uk) = ^^^(2^)'™^2„. (^^ ) , (0 < A: < 2n). (90) 



J«(fcAa;) - UO) ^ ^ l^ ^ ^\, B2rr.^2ikAxf. (91) 



Then 

4(2m-2)! 
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For example, the interaction (1801) with s = 3 gives 

Uk) - JM = -^ [k' - 4vrP + Arc'e] ^ -^ k'. (92) 

For s = 0, we have (Ref. [16], Sec. 5.4.2.9) the relation 

f;^^^^ = -ln[2sin(fc/2)]. (93) 

n=l 

Then, the limit Ax -^ gives 

JaikAx) ^ - ln(A;Ax) -^ oo. (94) 

For even numbers s, 

Ja{kAx) - 4(0) / \kAx\' -^ oo (95) 

since the expression has the logarithmic poles. 

The transition to the limit Ax — !► in Eq. (ISTj) with s = 1 gives 

^^^1^ - Gi A: u(A;, t) - J^{F {u{x, t))} = 0, (96) 

where Gi = ngAx is a finite parameter. The inverse Fourier transform of fl96|) leads to 
the continuous medium equation with coordinate derivative of first order: 

52 d 

T^u{x, t) - iGi — m(x, t)-F {u{x, t)) = 0. (97) 

at^ ox 

This equation can be considered as the nonlinear Schroedinger equation. 
The limit Ax ^ in Eq. flHTl) with s = 2m — 1 {m = 2,3, ...) gives 



d^u{k,t) 
where 



G2k^uik,t)-J^{F{u{x,t))} = 0, (98) 



\m— 1 /o-T^\2m— 2 



^ ( ^r (27rj , 

4(2m-2)! ^^ ^yy ! 

is a finite parameter. The inverse Fourier transform of (pHl) leads to the partial differential 
equation of second order: 



92 ^2 



u{x, t)-G2T^ u{x, t)-F {u{x, t)) = 0. (99) 



This equation can be considered as a nonlinear wave equation. 
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8 Noninteger power-law interaction 

Let us consider the power-law interaction with 

J(n) = \n\-^'-^'\ (100) 

where s is a positive noninteger number. 

Proposition 6. The power-law interaction U00\) with noninteger s is a-interaction 
with a = s for < s < 2, and a = 2 for s > 2. For < s < 2 (s ^ 1), the transform 
operation T maps the discrete equations with the interaction UU(J\) into the continuous 
medium equation with fractional Riesz derivatives of order a. For a > 2 (a ^ 3,4,5, ...), 
the continuous medium equation has the coordinate derivatives of second order. 

Proof. From Eq. fl^Ul) . we obtain the equation for u{k,t) in the form 

^^^^ + g lUkAx) - UO)] Hk, t) - T^{F K(t))} = 0, (101) 

where 



\n\ 

n=—oo ' ' 



For fractional positive a, the function (I102p can be presented by 

+ 00 -. 

UkAx) =J2^^ (e-*'^"^" + e*"^^) = Lt,+^{e''^n + I^^i+„(e-^'^^'), (103) 



n=l 



where Lip{z) is a polylogarithm function. Using the series representation of the polylog- 
arithm 1^: 



L.^(e-) = r(l - P){-zf-' + f; ^^^^, ""K ^ \z\<2n, /3^ 1,2,3..., (104) 



n=0 

we obtain 



UkAx) = A^ \Axr |fcr + 2 f; ^^^ + " ^""^ [Axf^-er, a y^ 0,1,2,3..., (105) 



n=0 

where ((z) is the Riemann zeta-f unction, |/i;Ax| < 27r, and 



A^ = 2 T{-a) cos f — ) . (106) 



From (11051) . we have 

Ja(0) = 2C(l + a). 
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Then 



Uk^^) - UO) = A^ \Axr \kr + 2 f2 ^^^^7^y!^(^^)'"(-^')"' (107) 



n=l 



where a 7^ 0, 1, 2, 3..., and |A;Ax| < 27r. 

Substitution of (ITOTI) into Eq. (ITOT]) gives 

n=l ^ '^ ' 

We will be interested in the limit Ax — > 0. Then Eq. (11081) can be written in a simple 
form 

— M(fc,t) + G, T,,A(A;)M(A;,t)-J-A{FK(t))} = 0, a ^0,1,2,..., (109) 

where we use the finite parameter 

G, = ^|Ax|™"^°'2>, (110) 

and 

. rA„|A;r-|Ax|2-"C(a-l)A;^ < « < 2, (« ^ 1); 

TaA[k] = < (111) 

[|Aa;|°-2A„|A:|°-C(a-l)fc2, a > 2, (a^3,4,...). 

The expression for 7^,a(^) can be considered as a Fourier transform of the interaction 
operator ([2]). From (II 101) . we see that g ^ 00 for the limit Ax -^ 0, and finite value of 

Ga- 

Note that (II 111) has a scale /cq: 

fco = |A,/C(« - l)|^/('~")|Ax|-^ (112) 

such that the nontrivial expression 7^,a(^) ~ |^|" appears only for < a < 2, (a 7^ 1), 
k < fco- 

The transition to the limit Ax -^ in Eq. (I109p gives 

-—u{k,t) + GMk)u{k,t)-J^-'{F{u{x,t))} = Q («^ 0,1,2,...), (113) 



where 



AJ*r 0<a<2. a^l; 

^ -C(«-1)P, 2<«, «^3,4,.... 



The inverse Fourier transform to (11131) is 

—^u{x,t) + GaZ,ix)u{x,t)-F{u{x,t)) = a ^0,1,2,..., (115) 

where 

\C(«-l)9V9|x|2, (a>2, «^3,4,...). 

Here, we have used the connection between the Riesz fractional derivative and its Fourier 
transform P: 

The properties of the Riesz derivative can be found in [H El [3l H] . 
As the result, we obtain the continuous medium equations 

Q2 Qa 

—u{x,t)-GaA^-^-—u{x,t) = F{u{x,t)), 0<a<2, (a ^ 1), (117) 

and 

—u{x,t) + G^C{c^-l)—-^u{x,t) = F{u{x,t)), a>2, (a ^3,4,...). (118) 
Analogously, the continuous limit for the system 

„ +00 

^ = 9Y1 \n-mr-'[un-Um]+F{un) (119) 

m=— oo 

gives the partial differential equations 

-u{x,t)-Go.A^——u{x,t) = F{u{x,t)), 0<a<2, [ay^l), (120) 

and 

— M(x,t) + G,C(a-l)^|^M(x,t) = F(M(x,t)), «>2, (a ^3,4,...). (121) 

For F{u) = 0, Eq. (I120p is the fractional kinetic equation that describes the fractional 
superdiffusion [331 [Ml [35]. If F{u) is a sum of linear and cubic terms, then Eq. (11201) has 
the form of the fractional Ginzburg-Landau equation [361 133 [Ml [39l HO]. A remarkable 
property of the dynamics described by the equation with fractional space derivatives is 
that the solutions have power-like tails. 
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9 Nonlinear long-range interaction 

In this section, we consider the discrete equations with nonhnear long-range interaction: 

+00 
in{u)= ^ J„(n, m)[/(u„) -/(Mm)], (122) 

m=— 00 

where f{u) is a nonhnear function of Un{t), and Ja{n,m) defines the a-interaction. As 
the example of Ja{n, m) = Ja{n — m), we can use the functions 

(-1)" 
'^'^^''^ " r{l + a/2 + n)r{l + a/2 -ny ^^^^^ 

We consider the interaction with f{u) = v^ and f{u) = u — gu^ that give the Burgers, 
Korteweg-de Vries and Boussinesq equations in the continuous limit for a = 1,2, 3, 4. If 
we use the fractional a in Eq. (I123p . we can obtain the fractional generalization of these 
equations. 

Proposition 7. The Fourier series transform JF^ maps the equations of motion 

-^ = 9 Yl Un-m)[fM-fiu,,)]+Fiu^), (124) 

m=— 00 

where F is an external on-site force, into the equation 
d^u{k,t) 



dt^ 



9[Ja{^) - J„(A;Ax)]^A{/(«n)} + ^A{i^(M„,)}, (125) 



where u{k,t) = J^A{un{t)} , and Ja{kAx) = J^A{J{n)}. 

If Ja{n) defines the a -interaction, then the continuous limit Ax -^ and the inverse 
Fourier transform give the equation 

d'^u(x t) d" 

^-^^ = G^A^^f{u{x, t)) + F{u{x, t)), (126) 

where Ga = g\Ax\'^ is a finite parameter. 

Proof. The Fourier series transform of the interaction term (11221) can be presented as 



-co +00 



J2 e-'="^^J„(n) = J2 E e-^^"^V(n,m)[/(«„) -/(«„)] 



n=— oo m=—oo 
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hoo +00 +00 +00 



= E E e-''="^V(n,m)/K)- J^ E e-''="^V(n,m)/K). (127) 

n=— 00 m=— 00 n=— 00 771=— oo 

For the first term on the r.h.s. of fll27p : 

+00 +00 +00 +00 

E E e-^^"^V(n,m)/K)= Y. ^-^'"^VK) $^ ^K) = ^a{/K)} J.(0), 

n=— 00 m=— 00 n=— 00 m,'=— 00 

(128) 
where we use J(?7i' + n, n) = J{m'). For the second term on the r.h.s. of 01271) : 



E E e-^'="^V(n,m)/(«™)= Y. /(^™) E e-^'='^'^V( 



n,m] 



-oo +00 



= Y /K)e"''™^" E e-''"'^V(n')=-FA{/(M„)} J„(A;Ax), (129) 

m=— 00 n'=— 00 

nVO 

where we use J{m, n' + m) = J{n'). 
As the result, we obtain Eq. (I125p . 
For the limit Ax — *> 0, Eq. (11251) can be written as 

—u{k, t) - G„ %,A{k) u{k, t) - T^{F {un{t))} = 0, (130) 

where we use finite parameter G^ = glAxl", and 

fa,A{k) = -^al^r - Ra{kAx)\Ax\-''. (131) 

Here, the function R^ satisfies the condition 

R^{kAx) ^ 
lim — — — I = 0. 

Ax-^O I Ax I" 

In the hmit Ax — i> 0, we get 

—u{k, t) - G„ %{k) T{f (m(x, t))} - T{F (m(x, t))] = 0, (132) 

where 

u{k,t) = Cu{k,t), fa{k) = Cfa,A{k) = -A„|A;|". 

The inverse Fourier transform of (11321) gives 

— m(x, t) - Ga T,(x) fiu{x, t)) - F {u{x, t)) = 0, (133) 
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where Ta{x) is an operator 



Ux) = J'-'iUk)} = A^^. (134) 



As the result, we obtain the continuous medium equation ([E 



Let us consider examples of quadratic-nonlinear long-range interactions. 

1) The continuous limit of the lattice equations 

-—— = gi 2_^ Ji{n,m)[u^-u^\+g2 2_^ J2{n,m)[Un - Um\, (135) 

m=—oo m=—OD 

where Ji{n) (i = 1, 2) define the aj-interactions with ai = 1 and 0:2 = 2, gives the Burgers 
equation [17] that is a nonlinear partial differential equation of second order: 

d d d'^ 

— M(a;, t) + ^1^(0;, t)—u{x, t) - G2-^u{x, t) = 0. (136) 

It is used in fluid dynamics as a simplified model for turbulence, boundary layer behav- 
ior, shock wave formation, and mass transport. If we consider J2{n,m) with fractional 
a2 = a, then we get the fractional Burgers equation that is suggested in Ref. |1H]. 

2) The continuous limit of the system of equations 

rj /,N +00 +00 

g =91 X^ Mn^™)[^l-Um]+93 ^ J3{n,m)[un-Um], (137) 

m=— oo m=—OD 

where Ji{n) (z = 1, 3) define the ttj-interactions with ai = 1 and a-^ = 3, gives Korteweg-de 
Vries equation 

d d d^ 

— m(x, t) - Giu{x, t)—u{x, t) + Ga— jm(x, t) = 0. (138) 

First formulated as part of an analysis of shallow-water waves in canals, it has subse- 
quently been found to be involved in a wide range of physics phenomena, especially those 
exhibiting shock waves, traveling waves, and solitons. Certain theoretical physics phe- 
nomena in the quantum mechanics domain are explained by means of a KdV model. It 
is used in fluid dynamics, aerodynamics, and continuum mechanics as a model for shock 
wave formation, solitons, turbulence, boundary layer behavior, and mass transport. 

If we use noninteger Oj-interactions for Jj(n), then we get the fractional generalization 
of the KdV equation 
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3) The continuous limit of 

d'^uJt) ^°° ^°° 

Z2 = 92 ^ Mn,m)[f{un) - fium)]+g4 5Z Mn,m)[un-Um], (139) 

m=— 00 m=— 00 

m^n rn^n 

where 

f(u) = u- gu"^, 

and Ji{n) define the aj-interactions with 0:2 = 2 and a^ = 4, gives the Boussinesq equation 
that is a nonhnear partial differential equation of fourth order 

f)2 f)2 f)2 f)A 

g^u{x, t) - G2g^u{x, t) + gG2^y{x, t) + G,—u{x, t) = 0. (140) 

This equation was formulated as part of an analysis of long waves in shallow water. It 
was subsequently applied to problems in the percolation of water in porous subsurface 
strata. It also crops up in the analysis of many other physical processes. 



10 Fractional derivatives from dispersion law 

Let us consider the three-dimensional lattice that is described by the equations of motion 

-^ = gYl -/(n,!!!) K-n„,] + F(Mn), (141) 

m=— oo 
rriT^n 

where n = (ni,n2, n^), and J(n, m) = J(n — m) = J(m — n). We suppose that Un{t) are 
Fourier coefficients of the function u(k,t): 

u{k,t)= Y, Wn(t)e-^'^'-" = J-AK(t)}, (142) 

n=— oo 

where k = {ki,k2, ks), and 



^. = E 



Here, aj are translational vectors of the lattice. The continuous medium model can be 
derived in the limit |aj| -^0. 

To derive the equation for u(k,t), we multiply ( ]14ip by exp(— ikrn), and summing 
over n. Then, we obtain 



du{k,t) 



J,(0)-J„(ka) M(k,t)+^A{F(Mn)}, (143) 
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where J^A{-^(wn)} is an operator notation for the Fourier series transform of F{u,^), and 

+00 
J,(ka) = Y, e"'''"" Jin). (144) 

n=— 00 

For the three-dimensional lattice, we define the a-interaction with a = (ai, 02, as), as 
an interaction that satisfies the conditions: 

lim j^"^'^,^ ~/°^°^^ = A.. (^ = 1>2,3), (145) 

where < |v4q,-| < 00. Conditions 0145p mean that 

3 3 

J,(0) - JM = J2^»M''' + 5^i?.,(k), (146) 

where 

lim R^X^)/\h\"^ = 0. (147) 

In the continuous limit (|aj| ^0), the a-interaction in the three-dimensional lattice gives 
the continuous medium equations with the derivatives d"'^ /dx"^ , d"'^ /dy"^ , and d"^ /dz°'^. 
Let us recall the appearance of the nonlinear parabolic equation [5ll [52l [53l ^^. 
Consider wave propagation in some media and present the wave vector k in the form 

k = ko + K = ko + K|| + K^, (148) 

where ko is the unperturbed wave vector and subscripts (||,-L) are taken respectively to 
the direction of ko. A symmetric dispersion law 

u{k) = u{k) = J„(ka) - J«(0) (149) 

for K = |k — ko| ^ fco = |ko| can be written as 

uik) = uji\k\) = uj{ko + [|k| - A;o]) ^ iuiko) + Vg (|k| - ko) + ]^v'g (|k| - ko)\ (150) 

where 

fduj\ , fd^Lu\ , , 

and 



|k| = |ko + k| = -y/(ko + fi;||)2 + fi;^ ^ /c^ + /tj| + —-k\. (152) 

ZKo 



Substitution of (IT52D into (IT50|) gives 
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where ujq = uj{kQ). Expressions (11431) and (11531) in the dual space ("momentum represen- 
tation") correspond to the following equation for u = u{r,t) in the coordinate space 

with respect to the field u = u(t,x,y,z), where x is along kg, and we use the operator 
correspondence between the dual space and usual space-time: 

.(k) ^ Z-, .,, ^ -.-, 

Equation (I154p is known as the nonlinear parabolic equation [211 [S21 ESI Ell • The change 
of variables from (t, x, y, z) to (t, x — Vgt, y, z) gives 

du V v' 

- i^ = ^^^« + f ^ll« - ^oM - F{u), (156) 

which is also known as the nonlinear Schrodinger equation. 

Wave propagation in oscillatory medium with long-range interaction of oscillators can 
be easily generalized by rewriting the dispersion law (11531) . in the following way: 



u 



{k) = u^ + VgK\\+G^{nlr'^ + Gp{i^l:{)^'\ (l<a,/3<2) (157) 



with new finite constants Ga-, and Gj3. 

Using the connection between Riesz fractional derivative and its Fourier transform [1] 

(_A^)-/2 ^^ {h.ITI\ (-A||)^/2 ^^ {nl;{)^'\ (158) 

we obtain from (11571) 



du . du 



I—— = —IV 



dt dx 



+ Ga(-Ax)"/'M + G^(-A||)^/2^ + ujou + F{u), (159) 



where u = u{t, x, y, z). By changing the variables from (t, a;, y, z) to (t, C,,y,z), C, = x — Vgt 
and using 

d\^~W\ 
we obtain from (11591) 



^\\f" = 4zr.=47r,^ (160) 



^^ = Go^i-A^r/\ + Gpi-A\\f/\ + uou + F{u), (161) 

which can be called the fractional nonlinear parabolic equation. For G^ = and F{u) = 
b\u\'^u, we get the fractional Ginzburg-Landau equation [36 | [371 [38 | l39l BO] . 
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We may consider one- dimensional simplifications of Eq. (11611) . i.e., 



where u = u(t, 0)^ = 2; 



du 




!gt, or 





(162) 



Oil d^ii 

^^ = G.^^+^o« + F(«), (163) 



where u = u(t,z). 

Let us comment on the physical structure of (11611) . The first and second terms on the 
right-hand side are related to wave propagation in oscillatory medium with long-range 
interaction of oscillators. The term with F{u) on the right-hand side of Eqs. (I159p . and 
(11611) correspond to wave interaction due to the nonlinear properties of the media. Thus, 
Eq. (I16ip can describe fractal processes of self-focusing and related issues. 

11 Conclusion 

One-dimensional system of long-range interacting oscillators serves as a model for numer- 
ous applications in physics, chemistry, biology, etc. Long-range interactions are important 
type of interactions for complex media. An interesting situation arises when we consider 
the wide class of a-interactions, where a is noninteger. A remarkable feature of these 
interactions is the existence of a transform operation that replaces the set of coupled indi- 
vidual oscillator equations by the continuous medium equation with the space derivative 
of noninteger order a. Such transform operation is an approximation that appears in the 
continuous limit. This limit allows us to consider different models in a unified way by 
applying tools of fractional calculus. 

Periodic space-localized oscillations, which arise in discrete systems, have been widely 
studied for short-range interactions. In the paper, the systems with long-range interac- 
tions were considered. The method to map the discrete equations of motion into the 
continuous fractional order differential equation is developed by the transform operation. 
It is known that the properties of a system with long-range interaction are very different 
from short-range one. The method of fractional calculus can be a new tool for the analysis 
of different lattice systems. 
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Appendix: Divergence of non-invatiant interaction term 

Noninvariant interaction term leads to the infinity in the continuous medium equation. 
To demonstrate this property, we prove the following proposition. 

Proposition 8. The a -interaction term 

+00 
g ^ J{n,m) Um, (164) 

m=— 00 
rriy^n 

where J{n,m) = \n — -ml"^""*"^) is not translation-invatiant. The transform operation T 
of the term (16^ leads to the divergence of order |Ax|~" in the continuous medium equa- 
tions. 

Let us prove this proposition for < a < 2 (a 7^ 1), and the following equations of 
motion 

-^+9 ^ J{n,m) u„i- F{un) = 0. (165) 

m=— 00 



Since 



hoc 



Y^ J{n,m)= Y^ |n-m|-("+^) ^0, 



m=—oo m=--oo 

then the interparticle interaction term in (11651) is noninvariant with respect to translations. 
To derive the equation for u{k,t), we multiply Eq. f ll65p by exp{—iknAx), and summing 
over n. Then, we obtain 

^^^ + g MkAx) u{k, t) - J^a{F {un{t))} = 0, (166) 

where Jaik) is defined by (I102p . Using (11051) . we present Eq. (I166p in the form 

^4§^ + g A^\Ax\" l^r u{k, t) + 2gC{a + l)u{k, t) + 
ot^ 

+ 2^ E ^^^^|^^yf^(Aa;)2"(-A;^)"u(A;, t) - ^a{F («„(t))} = 0, (167) 

where C, is the Riemann zeta- function and A^ is defined by (11060 . For the limit Ax -^ 
and < a < 2 (a 7^ 1), Eq. (I167P can be written as 

— n(A;, t) + G^ A,| A;]" u{k, t) + 2gC{a + l)u{k, t) - Ta{F (n„(t))} = 0, (168) 
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where < a < 2, a 7^ 1, and G^ = g\Ax\°' is a finite parameter. Note that g ^ 00 for 
Ax — >■ 0, if Ga is a finite. Therefore, the transition to the hmit Aa; ^ in Eq. (11681) 
gives the divergence term 

hm gC{a + l)u{k, t) = ({a + l)Gau{k, t) hm | Axj"" -^ 00. (169) 

To have the continuous model equations without divergences, we must consider [um{t) — Un{t)] 
instead of Um{t) in the interaction terms 
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